Constrained Molecular Dynamics Simulations of Atomic Ground-States 
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Constrained molecular dynamics(CoMD) model, previously introduced for nuclear dynamics, has 
been extended to the atomic structure and collision calculations. Quantum effects corresponding to 
the Pauli and Heisenberg principle are enforced by constraints, following the idea of the Lagrange 
multiplier method. Our calculations for small atomic system, H, He, Li, Be, F reproduce the ground- 
state binding energies reasonably, compared with the experimental data. We discuss also the shell 
splitting which is expected as a consequence of the e-e correlation. 



PACS numbers: 34.10.+x; 31.15.Qg 

Molecular dynamics approach is a powerful tool to sim- 
ulate nuclear heavy ion collisions and atomic collisions, 
due to its simplicity and ability to take into account the 
influence of correlations and fluctuations. However, as 
it is seen in the case of classical trajectory Monte Carlo 
method 1] , truly classical atoms, without constraints due 
to the Heisenberg and Pauli principle, are unstable. To 
describe the ground-state properties of the systems with 
molecular dynamics approach, the pseudo potential is 
often introduced to simulate the effects of the Heisen- 
berg uncertainty principle and the Pauli exclusion prin- 
ciple 0, ■ The method with pseudo potential is known 
as Fermionic molecular dynamics(FMD) and it is applied 
to the studies of various atomic processes, in which a fully 
quantum mechanical dynamical simulation encounters 
numerical problems difficult to overcome. To give some 
actual examples, the atom- ion collisions £|, the atomic 
ionization by laser fields Q , the capture of muons by hy- 
drogen atoms |(J , the formation of antiprotonic atom pfl 
and the ionization cross section calculations Q , all these 
processes have been investigated using molecular dynam- 
ics approach. In particular, Geyer and Rost Q proposed 
the quasi-classical calculation of the ionization cross sec- 
tion using classical propagation scheme. As the initial 
state of the target, they use the phase space distributions 
of the bound electrons which are directly obtained from 
quantum mechanical wave functions through the Wigner 
transformation. There they mention the importance of 
the energy spread of the phase space distributions. 

Meanwhile, a constrained molecular dynamics (CoMD) 
approach has been proposed to treat fermionic properties 
of the nucleons in nuclei . The approach has been suc- 
cessfully applied to study the Equation of State of the 
quark system as well [Toj . In this approach, the Pauli ex- 
clusion principle is accomplished by restricting the phase 
space occupancy /* to values less or equal to 1 [llj . The 
equation of motion with the constraints for each electron 
is derived on the basis of Lagrange multiplier method 
for constraints. The constraint of CoMD approach is 
thought as an alternative to the pseudo potential and 
can be easily extended to the case of the Heisenberg un- 
certainty principle as well. The constraints play the role 
of a "dissipative term" in the classical equation of motion 
and lead the system to its ground-state. But in CoMD, 



at variance with FMD approaches, the "dissipative term" 
can increase or decrease the energy of the system depend- 
ing on the phase space occupation. 

In this brief report, we apply CoMD to atomic systems 
for the purpose of determining their ground-states con- 
figurations. Particularly, we discuss some properties of 
ground-states atoms, i.e., binding energies (the total elec- 
tronic energies) and radial positions of the bound elec- 
trons. We discuss also the energy spread (variance of the 
binding energy) and radial variance of each bound elec- 
tron as well. In our approach, we prepare the ensemble 
of initial configurations. The binding energies and radial 
positions of the electrons in the systems are calculated 
as averaged values over the ensembles. Our approach is 
sufficient to obtain stable atomic ground-states provid- 
ing their atomic energies fairly accurately. Using the ob- 
tained ensembles of initial states which occupy different 
points in the phase space, molecular dynamics simula- 
tion with constraints for the atomic collision has been 
performed and applied to nuclear fusion enhancement 
factor calculations of D+rf reaction for astrophysical in- 
terests |12(. 

We describe the essence of CoMD briefly. In classi- 
cal molecular dynamics(CMD) one solves the Hamilton 
equations, i.e.,: 



dt 



= -V r t/(r,), 



(1) 



where we use relativistic kinematics: 



\Zpfc 2 + mfc i , TO, and 



U{vi) = 



N 

E • 



(2) 



are the energy, the mass and the potential of the i-th 
electron, respectively. Here, qi is the charge of the elec- 
tron i and Qj is the charge of electron(for 1 < j < N) 
or nucleus(for j = 0). The total Hamiltonian is written 
down as H(r; p) = (£i + £7(r) — mc 2 ). 

The shortcoming of the approach is the lack of the 
Pauli exclusion principle. CoMD gives a possibility to 
overcome the shortcoming. To take the feature of the 
Pauli blocking into account, we use the Lagrange mul- 
tiplier method for constraints. Our constraints which 



2 



correspond to the Pauli blocking is fi < 1 in terms 
of the occupation probability and can be directly re- 
lated to the distance of two particles, i.e., rijPij, in the 
phase space. Here = |rj — t*j | and pij = |p,: — Pj|. 
The relation fi < 1 is fulfilled, if r^pij 
where £p — 27r(3/47r) 2 / 3 2 1 / 3 , i,j refer only to electrons 
and Si,Sj(= ±1/2) are their spin projection. We can 
easily extend the approach to the Heisenberg principle 
where the constraint is expressed as TijPij > where 
£h = 1, i and j refer to the electrons and the nucleus. 
Using these constraints, the Lagrangian of the system 
can be written as 



A' 



£ = £p i .r i - J ff(r;p)+ £ A^ 



H I r ijPij 



- 1 



E x P ( r ijPij&s i ,s j \ 

n \~tp*~ ) 



(3) 



where Af and Xf are Lagrange multipliers for Pauli and 
Heisenberg principle respectively. The variational calcu- 
lus leads to: 



dr 1 = p j J 1 V fAf Af 
eft ft ft 4^ Uh £p 



ft 



op* 



(4) 



(5) 

From physical considerations we expect that the Pauli 
principle is stronger than the Heisenberg principle for the 
two closest identical electrons in the phase space, i.e., the 
particles i and for which r^Pij is smallest. While the 
Heisenberg principle must be enforced especially among 
the electrons and the nucleus. Thus we restrict the sum- 
mations in eqs. 10-© to those particles only. 

In order to obtain the atomic ground-state configura- 
tion, we perform the time integration of the eqs. Q and 
©. The values of Af and Af are determined depending 
on the magnitude of r^pij. If r^Pij is (smaller)larger 
than A has positive(negative) sign, changing the 

phase space occupancy of the system. The constraints 
work as the "dissipative term" in the case of the pseudo 
potential approach and lead automatically to the mini- 
mum energy, i.e., the ground-state of the system. The 
difference being that in the case of the model with the 
pseudo potential, a dissipative term decreases the total 
energy. In our case the total energy decreases or increases 
depending on the phase space occupancy. 

The electron configurations at the beginning of the 
time integration are prepared in the following way. In 
the case of even number of bound electrons, we locate 
a pair of them at the opposite points respect to the nu- 
cleus in the phase space. In this way, the center of mass 
of the electrons coincides with the position of the nucleus, 
i.e., total momentum of the electrons is zero. We do the 
same procedure for the odd-number electrons atom, ex- 
cluding an electron which is the outermost. Thus, at the 



beginning of the time integration we have an ensemble 
of electron configurations which occupy different points 
in the phase space microscopically. The integration of 
the eqs. @ and J5J) is performed using Hermite inte- 
gration scheme which is efficient and enables integration 
with high precision. The scheme adopts variable and in- 
dividual time-steps for each electron |l3j . Considering 
the nucleus rest frame, The binding energy of the atom 
is determined by 



^ nev 

B.E.= Vff(r;p), 



(0) 



where nev is the number of the events in the ensemble. 
The radial coordinates of the electron i (Ri) is given by 



H = — E vTr»: 
nev * — ' 



(7) 



The binding energies and the radial coordinates of the 
electron i of the atoms are calculated as an averaged value 
over the ensemble of events. 

We have applied the model to hydrogen, helium, 
lithium, beryllium, fluorine atoms. Fig. ^ and Fig. [2] 
show that the systems converge to their ground-state in 
the illustrative cases of lithium and beryllium atoms, re- 
spectively. The top panels show the time development 
of the average of rijPij/^h over all pair of particles and 
over events in the ensemble, we write it as ArAp/£?i. 
The middle panels and the bottom ones show the aver- 
age of binding energy and the radial coordinates of the 
electrons, respectively, over events. In the bottom panels 
each line corresponds to the radius of each bound elec- 
tron. Due to the constraints these values oscillate as a 
function of time and converge after some time. We de- 
termine the binding energy and the radius by taking the 
average over not only events but also over time. 

We summarize our results of the ground-state ener- 
gies and the radial coordinates of the bound electrons for 
small atomic systems in Table III and ITT1 respectively. For 
the purpose of utilizing the atomic configuration to colli- 
sion calculations, the comparison between our results and 
the ones from quantum mechanical Hartree-Fock(HF) is 
suggestive. In Table H] together with the ground-state 
energies from our method, results from the FMD [T3 |. 
HF [l^ methods and experimental values |l(ij | are shown. 
Since we determine the binding energies as an average 
of many events, our results have a variance \J (AB.E.) 2 , 
also included in the tables. We obtain ground-state bind- 
ing energies in good agreement with experimental data 
within variances. In Table ^ we compare our results 
of the radial coordinates of the electrons (Ri) for each 
atoms with FMD and HF methods. The results from 
FMD which are obtained using different parameter sets 
are shown in two columns(FMD and FMD 1 0). 
The column FMD 1 is with the optimized parameter 
sets. Note that our method gives smaller values as Ri 
than those from HF method. However one should no- 
tice that while comparison of our approach to FMD is 
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FIG. 1: The convergence to the ground state of the lithium 
atom. The values are averaged over 18 events. 



direct since they are all semi-classical molecular dynam- 
ics, thus the definition given in eqs. 0, (0) are exactly 
the same. On the other hand in the HF approach there 
are no "events" , the calculated values of HF are obtained 
through a smooth probability distribution, which is not 
the case for our model where we have (S-functions. Fur- 
thermore symmetries imposed to the system will be pre- 
served in HF, indeed they might be destroyed in CoMD 
because of correlations. We give an example of the Be 
case; initially we distributed the electrons pairwise, one 
opposite to the other respect to the nucleus both in the 
coordinate and momentum space, as stated above. This 
initial symmetry is broken in the simulation because of 
e-e correlation, nevertheless one can recognize two major 
shells within the variance ARi in the results. Certainly 
even in our calculations we can impose such a symmetry. 
In Tables [I] and ^ the numbers in parenthesis show in 
the case where we impose that electrons move pairwise 



opposite locations in the phase space, i.e., they are forced 
to keep the initial symmetries. Such a calculation gives a 
very similar B.E. as in the absence of forced symmetry, 
but now two major shells are clearly identified. From the 
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FIG. 2: The convergence to the ground state of the beryllium 
atom. The values are averaged over 20 events. 



comparison of these two cases it is obvious the splitting 
of the shells when relaxing the initial symmetry. 

We have presented results of Constrained molecular 
dynamics approach to describe the atomic ground-states 
configurations. We calculated the binding energies and 
the radial coordinates of the electrons in atoms. The total 
electronic energies for the ground-state atoms are given 
rather accurately, At last we stress that the intent of the 
CoMD simulation of the atomic systems is to applying it 
to the collision calculations and determining the Equa- 
tion of State of matter at very low temperatures where 
quantum effects play a decisive role. 

We thank Prof. J.S. Cohen for providing us his data. 
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TABLE I: Summary of the binding energies (in eV) of calculated systems. The binding energies with CoMD including error 
bars, with optimized FMD(FMD 1) with Hartree-Fock method |l5j| and experimental values [lf3 | are shown for each atom. 



CoMD 



FMD 1 



HF 



experimental 



H 

He 
Li 
Be 
F 



-13.3 ± 1.8 
-79.1 ± 6.9 
-196.8 ± 15. 

-371.2 ± 27. (-385.3 ± 29. 
-2560.3 ±192. 



-13.61 
-77.57 
-197.61 
-395.63 
-2408.3 



-13.61 
-77.87 
-202.26 
-396.56 
-2705.1 



-13.6 
-79.1 
-203.6 
-399.4 
-2717.4 



"with forced symmetry 



TABLE II: Summary of the radial coordinates of bound electrons(rc) (in atomic unit) of calculated systems. The rc with 
CoMD, with FMD(FMD |3 and optimized FMD 1 [3) and with Hartree-Fock method are shown. 

CoMD 



H 

He 

Li 

Be 



1.0 ± 0.2 

0.55 ± 0.09 
0.60 ± 0.1 

0.34 ± 0.1 
0.5 ± 0.2 
2.0 ± 0.6 
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(0.36 ± 
(0.36 ± 
(1.60 ± 
(1.60 ± 



0.17) 6 
0.17) 
0.34) 
0.34) 



FMD 


FMD 1 


HF 


1.0 


1.0 


1.0 (Is) 


0.5714 


0.6139 


0.9273(ls) 


0.5714 


0.6139 




0.3506 


0.3818 


0.5731(ls) 


0.3673 


0.3885 


3.8737(2s) 


1.4419 


4.0604 




0.2565 


0.2687 


0.4150(ls) 


0.2565 


0.2687 


2.6494(2s) 


0.9458 


2.9738 




0.9458 


2.9738 






0.2625 


0.1757(ls) 




0.2624 


1.0011(2s) 




0.2622 


1.0848(2p) 




0.2622 






0.2668 






0.2671 






0.6931 






0.6952 






0.7818 





'same as in Table |I]with forced symmetry 
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